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We study the ground state properties and the excitation spectrum of bosons which, in addition 
to a short-range repulsive two body potential, interact through the exchange of some dispersionless 
bosonic modes. The latter induces a time dependent (retarded) boson-boson interaction which 
is attractive in the static limit. Moreover the coupling with dispersionless modes introduces a 
reference frame for the moving boson system and hence breaks the Galilean invariance of this 
system. The ground state of such a system is depleted linearly in the boson density due to the zero 
point fluctuations driven by the retarded part of the interaction. Both quasiparticle (microscopic) 
and compressional (macroscopic) sound velocities of the system are studied. The microscopic sound 
velocity is calculated up the second order in the effective two body interaction in a perturbative 
treatment, similar to that of Beliaev for the dilute weakly interacting Bose gas. The hydrodynamic 
equations are used to obtain the macroscopic sound velocity. We show that these velocities are 
identical within our perturbative approach. We present analytical results for them in terms of two 
dimensional parameters - an effective interaction strength and an adiabaticity parameter - which 
characterize the system. We find that due the presence of several competing effects, which determine 
the speed of the sound of the system, three qualitatively different regimes can be in principle realized 
in the parameter space and discuss them on physical grounds. 
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The superfluid properties of Bose-systems have beep, 
studied and greatly understood since the early fifties n 
While an ideal Bose gas condenses but does not exhibit 
superfluidity, an interacting Bose gas with repulsive in- 
teraction does. 

The microscopic description of excitation spectrum of 
weakly interacting dilute Bose n gas (WIDBG) has been 
first considered by Bogoliubov.B ft has been shown that 
the interaction among the particles modifies the charac- 
ter of low lying elementary excitation in the condensed 
state of the systema and as a result long-wavelength ex- 
citations obey the linear dispersion law, thus having a 
sound-wave-like character rather then particle-like as in 
the normal state of the system. Within this lowest or- 
der treatment, referred to nowadays as Bogoliubov ap- 
proximation, the microscopic sound velocity was found 
to be identical to the macroscopic, compressional, sound 
velocity of the system. Going beyond this lowest order 
theory, Beliaev_has developed a perturbation expansion 
of the WIDBGEjEI and considered the second order cor- 
rections to the Bogoliubov results. His results confirmed 
the equivalence of microscopic and macroscopic sound 
velocities BlI This equivalence has been generalized to all 
orders of perturbation theory by Gavoret and Nozieres in 
Ref. |7] invoking explicitly the Galilean invariance of the 
system. 

In a recent paperi we have addressed ourselves to 
a class of systems composed of two coupled bosonic 
subsystems: propagating bosons, interacting through a 
short-range repulsive interaction, and local dispersionless 
bosonic degree's of freedom (for brevity hereafter termed 



as phonons ) to which they are coupled. The coupling 
to the latter not only supplements the short range two 
body potential of bosons by a time dependent counter- 
part (attractive at low energies) but also introduces a 
reference frame for the moving boson system and hence 
breaks the Galilean invariance of this system. The aim of 
the present work is to study the ground state properties 
and the excitation spectrum of such a system, focusing 
on the renormalization of long- wavelength sound-like ex- 
citation due to the presence of the competing interactions 
in the system mentioned above. 

We first determine the microscopic sound velocity from 
the single-particle excitation spectrum. The latter has 
been calculated in Ref. |^ within the Beliaev-PopovB'cl 
second-order perturbation theory generalized to our sys- 
tem. Here we consider the continuum version of the 
model discussed in Ref. || and show that the leading or- 
der terms in a density expansion can be treated analyt- 
ically. We discuss the modification of the perturbation 
expansion due to the presence of retarded interaction and 
the range of validity of our treatment. To calculate the 
macroscopic sound velocity we base ourself on the hy- 
drodynamic equations, explicitly incorporating the non- 
Galilean invariance of the system. We find that both 
velocities are identical up to the order considered in the 
present approximation. 

The Hamiltonian which describes the above discussed 
scenario is given by 

q k,k',p 
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+ UJ ^( a q a 1 + 2) ~ 51 6 k 6 k+q[aq + 0-q] , (1) 

q k,q 

where b^P and a q ^ denote the boson and phonon anni- 
hilation (creation) operators, e q = g 2 /(2m), m and /i 
is the boson mass and chemical potential, respectively. 
The interaction among bosons is assumed to be properly 
renormalized and is given by g = Aira/m. We note, that 
working with a renormalized interaction, one has to be 
careful to avoid double counting of the diagrams which 
have already been included in the ladder sum for the 
renormalized interaction, that defines the s-wave scat- 
tering length a of two particles in vacuum. The strength 
of the boson-phonon coupling is given by aujQ, (a being 
dimensionless parameter) and loq denotes the frequency 
of the dispersionless phonon mode. 

In our treatment we assume both, the dilute limit of 
the boson system and weak boson-phonon coupling in 
order to set up a perturbation theory for the above in- 
troduced Hamiltonian at T = 0. The phonon mode 
can be effectively integrated out leading to the effec- 
tive two body potential r q , w = g + a 2 uj 2 Vq}j, where 

f q °L = 2(jJo/(^ 2 — Wq) is the bare phonon Green's func- 
tion. The phonon meditated interaction vanishes at high 
frequencies as 1/uj 2 and hence no special renormaliza- 
tion or ultraviolet cut-off is required in order to perform 
calculations to a given order. 

The lowest order approach to the above Hamiltonian 
is constructed within the Bogoliubov schemeEl: by sep- 
arating out in Eq. ([!]) the condensate part of the Bose 

field 6k = frk + n>c $k,o (n c being the condensate fraction) 
and keeping only the bilinear part of the resulting Hamil- 
tonian. To this order of approximation one obtains two 
branches of the excitation spectrum of the systemp: 

<4 >q = ± ^l? + 16«2^ £q n c (2) 

with Eq = [e q + 2gn c £q] s denoting the spectrum of Bo- 
goliubov quasi-particles in the absence of the coupling 
with phonons. 

When the bosons condense, the phonons, being ini- 
tially coupled to the boson density (symmetry restoring 
variable in the case of gauge symmetry breaking), get 
hybridized with the Goldstone mode. The resulting two 
normal modes (^) describe Bogoliubov type excitations 
and Einstein phonons. For momenta close to where the 
level crossing of bare excitation spectra occurs neither 
mode is predominantly a phonon nor a Bogoliubov quasi- 
particle. In the long-wavelength limit the lower branch 
reduces to a sound- like dispersion w_. q ~ voq, where 
v o = [(.9 — 2a 2 uJo)n c /m} a is a characteristic sound veloc- 
ity with the repulsive interaction having being reduced 
by the static part of the attractive phonon mediated in- 
teraction. 

The Bogoliubov scheme amounts to taking into ac- 
count collisions between condensate quasi-particles, be- 
tween condensate and out-of-condensate quasi-particles, 



while neglecting scattering among out-of-condensate 
quasi-particles. The interaction among the particles de- 
pletes the ground state condensate leading to a finite 
density n of non condensed particles. For the weakly 
interacting dilute Bose gas, the depletion is given by 
n = n — n c ~ n 3 j3,,and is small in the low density limit 
n <C 1. BeliaevOla has considered the second ordered 
corrections to the Bogoliubov result and has shown that 
the expansion parameter of the theory is the so called 
gas parameter an 1 / 3 . The latter measures the effective 
range of the interaction a (the s-wave scattering length) 
relative to the inter-particle distance d ~ n 1 / 3 . All the 
second-order corrections to the Bogoliubov results (for 
example corrections to the chemical potential /i, squared 
sound velocity v 2 ) appear with higher powers in density 
(~ n 3 / 2 ) and are small in the dilute limit even if the 
interaction is strong. As it follows, these arguments no 
longer apply to the case we consider. The reason is the 
following: even in the vanishing density limit the system 
does not reduce to the non interacting one and one ends 
up with the well known polaron problem. In other words, 
due to the retarded character of phonon mediated inter- 
action, one particle can interact with itself, by creating 
the lattice deformation at given time and absorbing it 
later on. 

A generalization of the Beliaev-PopovoD second-order 
perturbation theory to the presently studied problem has 
been discussed by us in Ref. [| As we have already men- 
tioned, here we restrict ourselves to the leading order in 
density contributions, which allows us to obtain analyti- 
cal results. We calculate the corrections to the chemical 
potential, microscopic sound velocity and discuss the ex- 
pansion parameter and the limit of the validity of the 
perturbation theory. 




FIG. 1. a) Leading order in density contribution to ground 
state depletion, b) normal component of the boson self-energy 
due to the coupling with phonons, c) anomalous component 
of the boson self-energy, d) renormalized four point vertex. 



To begin with we calculate the depletion of the ground 
state h. It is given by the number of out-of-condensed 
particles and vanishes for the ideal Bose gas. In Fig. la 
we express graphically the linear in density contribution 
to the depletion. The straight and wave lines stand for 
the boson and phonon bare Green's functions, respec- 
tively, the dot represents the boson-phonon vertex and 
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the zig-zag line the condensate (given by a factor n]/ 2 ). 
This contribution is exclusively due to the phonon me- 
diated retarded interaction. It vanishes for the instanta- 
neous interaction since in this case the end points (see 
Fig. la) would correspond to the same time, which in- 
volves both "particle" and "hole" type excitations. Only 
the former one is non vanishing for T = 0. 

The corresponding analytical expression of the dia- 
gram shown in Fig. la is as follows 
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In Eq. 



D = 



(uj — e q + irj) 1 denotes the bear bo- 



son Green's function. As a result, to the leading order in 
density, the ground state depletion and the condensate 
fraction takes the form: 



7 A 
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1 + 7A/2 
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where we have introduced two dimcnsionlcss parameters: 
A = 2a 2 u>o/g which is the strength of phonon mediated 
interaction relative to the direct boson-boson repulsion 
and 7 = a/a,Q being the ratio of two characteristic lengths 
- the s-wave scattering length a and o,q = 1/qo with 
<7o = \Z2mujo being the momentum at which level cross- 
ing of the boson and the phonon bear spectrum occurs. 
We note that the parameter 7 always appear in the com- 
bination 7 A which does not explicitly involve the boson- 
boson repulsion and is equivalent to the so-called Migdal 
parameter in the electron-phonon problem. However, as 
it follows, the interpretation of our results in terms of 7 
and A is more transparent and we keep this notations, 
rather using the conventional one. 

The renormalization of the spectrum given by Eq. (||), 
due to the second order corrections can be obtained from 
the poles of the dressed boson propagator. The diag- 
onal (Gq,u) and off-diagonal (Gq u>) elements of the bo- 
son Green's function are expressed through the normal 
(Eq jU ) and anomalous (£q, w ) components of the self- 
energy 



Gq,u, = [W + £q + 



]/Eq ;[i ,, Gt\,u [ £q,aj]/Eq ?ti ;, 

], (5) 



where e q = e q - fi, Sq tU = [E qjU , + S q ,_ w - 2£ qi J/2 and 
Aq,u> — [£q i(il — £ q ._ w ]/2 are symmetric and antisymmet- 
ric functions in-w, respectively. Invoking the Hugenholtz- 
Pines theoremtl we have the relation fi = £0,0 — ^0.0 = 
So,o- 

Expanding the self-energies in the long- wavelength, 
low frequency limit as follows 



we obtain the microscopic sound velocity from Eq. (|^) as 

2 £o : o(l + '-'2) 



77l[(l + |^|) 2 + 2|5 1 |S , ] 



(7) 



Let us first discuss the self-energy expansion parame- 
ters defined in Eq. (^|). As we have already mentioned, 
we are interested in the leading order (linear in density) 
contributions. Since the anomalous self-energy £0,0 en- 
tering into Eq. (Q) vanishes linearly with density, there 
will be no contribution to this order coming from Si and 
only those from the remaining \A\ and <§2 (which has 
nonzero value in the limit of vanishing density) will con- 
tribute to this order. The first non-zero contribution to 
A comes from the diagram presented in Fig. lb. Substi- 
tuting the bare boson and phonon Green's functions in 
the corresponding formula of this diagrams, and expand- 
ing the result of the integration in powers of the external 
frequency u>, after straightforward algebra one obtains 
A = 7A/2. One also verifies, that for the dispersionless 
Einstein phonons considered here, the phonon mediated 
self-energy (see Fig. lb) being local in space is momen- 
tum independent. Hence, for the momentum expansion 
parameter 1S2 @ one has 1S2 = within the same order 
of approximation. 

At this point, it is interesting to check whether the 
exact relationsB 



|^| = dh/dn c , S 2 q 2 



(8) 



Aq :UJ = Au, Sq^ — 11 = S\bj 2 + S 2 q 2 



(6) 



derived in Ref. still holds, given the above obtained 
perturbative results. In Eq. (||) e' and e' q stand for the 
energy per unit volume of the system at rest and the 
system moving uniformly with a speed v = q/m, respec- 
tively (prime denotes that the contribution from the con- 
densate phase is eliminated). With the help of Eq. (Q), 
one can easily verify that the first relation Eq. (||) is 
recovered and not violated by the perturbative theory. 
As for the second relation, invoking the Galilean invari- 
ance of the system it has been linked in Ref. ^ to the 
change in kinetic energy of the moving system leading to 
62 = n~ 1 n/(2m). For the present case [see Eq. (|])] this 
gives a finite value 1S2 = 7 A/ (4m) instead of zero as has 
been discussed above. However, as it will be shown in 
details below, due to the non-Galilean invariance in the 
present case there is an additional contribution to the en- 
ergy of the moving system, stemming from the retarded 
boson phonon interaction. Considering both these con- 
tributions, it turns out that to within this order of ap- 
proximation of our perturbative treatment they cancel 
each other and one recovers 1S2 = 0. 

To the leading order in a density expansion the anoma- 
lous self-energy is given by the diagram shown in Fig. lc, 
corresponding to £0.0 = ^ n c with Y denoting a four point 
vertex dressed by the phonon mediated interaction (see 
Fig. Id) with all external momenta and frequencies taken 
to be zero. The renormalization of the four point vertex 



3 



in the second order is shown in Fig. Id. The diagrams d% 
and c?4 from Fig. Id represent the renormalization of the 
boson-boson and the boson-phonon vertices, respectively. 
The diagrams d§ and describe the two-body i-matrix 
renormalization arising from the phonon-mediated inter- 
action. The last two diagrams (^7 and d&) represent 
the exchange processes and, like the diagrams for vertex 
renormalization, are nonzero due to the time dependence 
of the phonon mediated interaction. We note that there 
are no equivalent diagrams with both interaction lines 
being the boson-boson renormalized interaction g, since 
either these contributions involve backward propagators 
of the bosons in the vacuum and hence vanish or they 
have already been included in g and do not have to be 
double counted. After substituting the bare boson and 
phonon Green's functions in the corresponding standard 
analytical formulas of the above discussed diagrams, all 
the integration can be done analytically and one arrives 
to the following expression for f (see Ref. |l(|): 



.g[l-A + 4 7 A-177A 2 /8] 



(9) 



Based on the above obtained result the sound velocity 
from Eq. (JTI) can be written as: 



gn 
m 



1 + -(-8 + 20 7 - 5 7 A) + <3(7 2 A 2 ) 
o 



(10) 



This expression implies the perturbation expansion in 
powers of A and 7A. The dilute limit does not necessary 
imply the convergence of the perturbation expansion and 
we have to assume moreover A << 1 and 7 << 1/A. 

Before providing a physical interpretation of these re- 
sults, let us examine the macroscopic sound velocity of 
the system to within the same order of perturbation the- 
ory. For that purpose we start with the well known equa- 
tions of motion for the two conjugate dynamical vari- 
ables: the number operator n and the macroscopic phase 
of the condensate 6. It reads as: 



On 



de 



n,(f> 



50 



de 



n,cj> 



dn 



(11) 



where e n ^ denotes the energy of the system per unit vol- 
ume. The next step is to consider a hydrodynamic poten- 
tial flow of the boson system with the superfluid velocity 
defined as v s = V</>/m. Then the change in energy of the 
system at finite but small v a can be written, without loss 
of generality, as 6e ntVs — e niVs — e n fl = Amv 2 /2, where 
A is the phase stiffness of the superfluid system. Taking 
the grad from both sides of the second equation (|Tl| ) and 
using the thermodynamic relation de n ^/dn — a one ar- 
rives to the Josephson relation for the superfluid velocity 
field mdv s /dt = — V/i. The two equations of motion ( p"l| ) 
can be combined into: 



d 2 n 



c 2 V 2 n, 



(12) 



with c 2 = A/x being the speed of macroscopic sound, 
and 1/x = d[ijdn the compressibility of the system. 



In the case of Galilean invariant system the change of 
the system's energy is exclusively due to the kinetic en- 
ergy term in the Hamiltonian and is given by 5e„ tVs = 
^ e n"u = nmv 2 /2. The potential energy, being only a 
function of the relative positions of the particles, does 
not change when the boson system moves with finite ve- 
locity. Then the phase stiffness is given by A = n/m 
and one recovers the well known expression for the speed 
of compressional sound c 2 = n/(mx)- However, in the 
case of coupling with phonons, the system is no longer 
Galilean invariant since then this coupling introduces the 
reference frame for the moving system of bosons. In other 
words, the phonon mediated interaction, being retarded, 
depends on the frequency and hence on the velocity of the 
particles. Therefore, there will be an additional contribu- 
tion to 8e n ,v s due to the boson-phonon coupling 5e^~ p . 




FIG. 2. Ground state energy correction due to the bo- 
son-phonon coupling. The quantities in the parentheses cor- 
respond to the frame moving together with the boson system. 

In Fig. 2 we present graphical expression for the first 
nonzero contribution to <5e^~ p . In the frame at which 
the phonons are at rest, the condensate mode has a fi- 
nite momentum equal to v s /m. To calculate <5e^~ p it is 
convenient to perform a Galilean transformation to the 
frame in which the boson system is at rest. This amounts 
to replacing the frequency u> appearing in the phonon 
propagator in the diagram shown in Fig. 2 by (lo — v s q). 
Evaluating this contribution and expanding the result 
for small v s one finds 5e®~ F = — jXnmv 2 / [2(2 + 7A)]. 
Combining the two contributions, 6e^™ s and <5e^~ p , we 
arrive at 6e n<V3 = A*m 2 v 2 /2, where the renormalized 
phase stiffness is given by 



A* =n/[(m(l + 7 A/2)]. 



(13) 



The next step is calculate compressibility of the system. 
To this end we evaluate the chemical potential by using 
the Hugenholtz-Pines relation u = £0.0 — ^0.0 B Examin- 
ing the various diagrams for Eo,o & n d £0,0 one finds that 
most of them cancel pairwise, and to leading order (linear 
in density) we have \i = const + (g — 2a 2 wo)n + /i^ 2 ) . The 
constant term describes the negative density independent 
shift of the chemical potential due to the coupling with 
phonons and does not contribute to the compressibility. 
The second term is the first order contribution given by 
the reduced boson repulsion arising from the static at- 
tractive part of phonon mediated interaction. The last 
term describes the second order corrections and can be 
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written as = [2d 5 + d 6 + 2dj + ds]n c (see Ref. [ll]) 
where di denotes the diagrams shown in Fig. Id and the 
corresponding values are given in Ref. [h]. The final result 
for compressibility is written as 



X 



1 = 1 - X - 37A - 9 7 A 2 / 



(14) 



Finally, upon using Eqs.@, @, and @ one ob- 
tains the macroscopic sound velocity: 



gn 
m 



X 



1 + -(-8 + 20 7 - 5 7 A) + 0(7 2 A 2 ) 



(15) 



which is thus identical to the microscopic one (|i0|). 

Let us now discuss our main findings. Our pertur- 
bative expansion assumes A C 1 and 7 A -C 1. The 
first condition is necessary to insure the stability of the 
condensate against the collapse when the attractive in- 
teraction overcompensates the repulsive one, the second 
is equivalent to the one of the Migdal approximation 
in the electron-phonon problem. For small A -C 1 one 



can neglect the last, higher order, term in Eq. (15) and 
obtains the following expression for the sound velocity 
c 2 ~ c 2 [l + A(— 1 + 57/2)], with Cq = gn/m being the 
speed the sound in the absence of any boson-phonon cou- 
pling. Analyzing this expression as a function of 7 one 
finds three different regimes of behavior: 

i) for 7 < 2/5, corresponding to a small phonon mode 
frequency, the sound velocity decreases when the boson- 
phonon coupling is switched on, 

ii) for 7 = 2/5 the linear in A contribution is exactly zero, 
and the sound velocity does not depend on A, provided 
the latter is small as we have already assumed, and 

iii) if 7 > 2/5, corresponding to a high phonon mode 
frequency, the coefficient in front of A becomes positive 
and the sound becomes stiffer when the boson-phonon 
coupling is switched on. 

While intuitively the softening of the sound would be 
naturally expected, the increase of the sound is a rather 
surprising result. It can be understood as follows. The 
softening of the sound is due to two effects: First the 
phonon mediated interaction, being attractive, reduces 
the boson-boson repulsion, and second it enhances the 
boson mass. However, if the reduction of the repulsive 
interaction in the lowest order is not strong enough, it 
could be overcompensated by the enhancing of the boson 
density vertex due to the coupling with phonons. More- 
over, this enhancement could also overcompensate the 
reduction of sound velocity due to the enhancement of 
the boson mass. Both effects, mass enhancement and 
vertex renormalization have the same origin, being only 
driven by the time dependence of the phonon mediated 
interaction. Hence these effects have the same order of 
magnitude and, acting in different directions, strongly 
compete. This is manifest in the three different regimes 
found in our perturbative approach. We would like to 
point out, that all regimes can in principle be realized 
within the range of validity of the present treatment. 



In our previous paper ,cl analyzing numerically Beliaev- 
Popov self energies for different values of a (boson- 
phonon coupling strength), without restricting ourself to 
the leading order terms in density expansion, discussed 
here, we came to the conclusion that the sound veloc- 
ity was practically unaffected by the coupling with the 
phonons. This situation which corresponds to the sec- 
ond regime discussed above. 

In conclusion, we have studied the ground state proper- 
ties and the excitation spectrum of the system of bosons 
which, in addition to short-range repulsive two body po- 
tential, interact through the exchange of the other disper- 
sionless bosonic modes, for brevity termed as phonons. 
We have assumed the dilute limit of the boson system 
and weak boson-phonon coupling and considered the sec- 
ond order perturbation theory, equivalent of the Beliaev- 
Popov treatment of the weakly interacting dilute Bose 
gas. We have shown that the ground state of such a 
system is depleted linearly in the boson density due to 
the zero point fluctuations driven by retarded part of the 
interaction. We have derived analytical expressions for 
the microscopic and macroscopic sound velocities of the 
system and shown that up to the order of our pertur- 
bative approach they are identical. We found that three 
physically different regimes can be realized in the param- 
eter space of the system and discussed them on physical 
grounds. While one can think of a number of phenomena 
which might be described in terms of the model we con- 
sidered, we leave the application of our study to realistic 
systems for some future work. 
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deriving these equations no Galilean invariance has been 
invoked in Ref. M and hence they still hold in present case. 
First we point out, that the diagram given in Fig. la 
that has been already evaluated Eq. (^) is the basic el- 
ement of the diagrams d$, (I4 and dy shown in Fig. Id. 
Based on this observation one finds d$ = dj — gyX/2 
and di = — 2q 2 ujo7A/2. One farther verifies, that dia- 
gram ds, which can be obtained from dj by reversing 
an arrow of one of the boson propagators, is twice big- 
ger. This leads to ds = 2dj = (77A since in the former 
one (ds) both two poles (in lower and upper half of the 
complex plane) of the phonon propagator contribute. The 
same holds for d§ and ds - and it is straightforward to 
show da = 2dg = — 3a 2 ajo7A/2. Putting all these contri- 
butions together and noting that the diagrams d^, di, ds, 
and dr enters with the symmetric factor 2, one obtains 



f = di + d2 + 8ds + 2c?4 + 3d$. Substituting in the last 
expression those corresponding values of the diagrams, dis- 
cussed above, one recovers Eq. (|^). 
11 The second order corrections to the chemical potential fi^ 
have been discussed in details in Ref. pi The diagrammatic 
representation of the contribution to jr 2 ' due to the boson 
phonon coupling is shown in Fig. 7 of that paper. The con- 
tributions due to the boson-boson interaction is obtained 
by replacing the phonon propagator by the boson-boson in- 
teraction line in these diagrams. The leading order in den- 
sity terms, can then be obtained by replacing the dressed 
boson propagators in these diagrams by corresponding per- 
turbation series and keeping only the leading, linear in den- 
sity terms. In this way the same diagrams as those shown in 
Fig. Id of the present paper are generated and one recovers 
the analytical formula for /i' 2 ' presented in the text. 



6 



